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We report on the observation of multiphoton interband absorption processes for quantum gases in 
shaken light crystals. Periodic inertial forcing, induced by a spatial motion of the lattice potential, 
drives multiphoton interband excitations of up to the ninth order. The occurrence of such excitation 
features is systematically investigated with respect to the potential depth and the driving amplitude. 

Ab initio calculations of resonance positions as well as numerical evaluation of their strengths exhibit 
good agreement with experimental data. In addition our findings could make it possible to reach 
novel phases of quantum matter by tailoring appropriate driving schemes. 
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I. INTRODUCTION 

Periodic driving of quantum systems allows for the tar¬ 
geted engineering of exotic properties. In recent years, 
this approach has been very successfully utilized in vari¬ 
ous fields of physics. While time-periodic forcing of solid- 
state materials yields access to, e.g., topological band 
structures [1-4] and high-Tc superconductors [5-8] it is 
also applied to trapped ions [9], photonic crystals [10] 
and for ultracold atomic ensembles [11, 12]. Quantum 
gases in optical lattices are particularly well suited as 
they are almost perfectly isolated from their environment 
and allow for unprecedented control in a time-dependent 
fashion. So far, experimental studies have focused on the 
creation of tunable artificial gauge potentials and large 
effective magnetic fluxes in optical lattices [13-22] , which 
allow for the observation of exotic phenomena such as ge¬ 
ometrical frustration [13], chiral Meissner currents [20], 
Ising magnetism [17, 18, 23], and topological band struc¬ 
tures [21, 22]. Further driving schemes have been pro¬ 
posed that aim for the realization of non-Abelian gauge 
fields [24, 25]. 

According to the Floquet theorem, the evolution of 
time-periodic systems can be described in terms of a 
time-periodic unitary operator and a time-independent 
effective Hamiltonian [26-29]. The underlying princi¬ 
ple of all driving schemes is that the properties of the 
driven system are determined by the effective Hamil¬ 
tonian, which might exhibit interesting novel features. 
This so-called Floquet engineering typically assumes that 
excited Bloch bands can be neglected. However, quan¬ 
tum gases in periodically driven optical lattices exhibit 
close analogies with laser-irradiated solid-state materials 
in which nonlinear processes play a crucial role at large 
field strengths [30-32]. Indeed, similar to an oscillat¬ 
ing light field, external periodic forcing of optical lattice 
systems is expected to induce significant multiphoton ex¬ 
citations between energy bands [33]. Thus, a deeper un¬ 
derstanding of such excitations processes is essential for 
tailoring appropriate driving schemes. 

Beyond its relevance to Floquet engineering, peri¬ 
odic forcing of optical lattices allows for the investiga¬ 


tion of multiphoton absorption (MPA) processes in well- 
controlled model systems, in which the band structure 
and the interaction strength are fully tunable. The be¬ 
havior of such non-linear excitations in interacting sys¬ 
tems opens far-reaching lines of investigation. 

In this paper we present a systematic study of mul¬ 
tiphoton excitations of ultracold bosonic quantum gases 
in driven optical lattices. Periodic inertial forcing of the 
atomic ensemble that is induced by a rapid shaking of the 
rigid lattice potential results in interband excitations due 
to MPA of low-energy driving photons. The emergence of 
resonance features corresponding to MPA between differ¬ 
ent Bloch bands is investigated with respect to the lattice 
depth and the driving field intensity in a one dimensional 
optical lattice. We extend our experimental studies to a 
two-dimensional triangular lattice in the regime of neg¬ 
ative effective tunneling. The positions of the observed 
resonances as well as their relative strengths are in good 
agreement with theoretical simulations. 

II. EXPERIMENTAL SETUP 

Here, multiphoton interband excitations are investi¬ 
gated in an ensemble of ultracold bosonic ^^Rb atoms 
that is confined in a red-detuned one-dimensional op¬ 
tical lattice. As depicted in Fig. 1(a), the lattice con¬ 
sists of a pair of running-wave laser beams of wavelength 
Xl = 830 nm that are arranged at an angle of 117.1° 
with respect to each other in the a;^-plane. The result¬ 
ing standing light wave Idat(^) = Vq cos(27rx/a)/2 has 
a lattice spacing of a = 486.5 nm. Its potential depth 
Vq is conveniently expressed in units of the recoil energy 
Free = ^^^i/(2M), denoting the kinetic energy that is 
transferred to an atom of mass M by absorbing a lat¬ 
tice photon with wavenumber kp = ‘I'KjXp, where h 
is the reduced Planck constant. In the following, the 
band structure of the lattice is denoted by ea{q)^ with 
band index a = 0,1,2,... and quasimomentum wave 
number q G {—nla^TTla] in the x-direction. The lowest 
bands are dominated by tunneling between neighboring 
lattice minima with tunneling parameters Jo > 0 and 
Ji < 0, respectively, so that So{q) ~ —2 Jq cos(ag') and 
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Lattice depth Vg (units of E^ec) ^ (units of n/a) 


FIG. 1. (Color online) (a) Setup of the driven ID lattice with 
both beams being linearly polarized along the z-axis. (b) 
Illustration of two-photon transitions into the first excited 
band at the minima at g = 0 for Jq^ > 0 (black) and q = 
TV/a for Jq^ < 0 (gray), (c) Multiphoton transition energies 
(white lines) to the first excited band according to Eq. (6) for 
K = 3.82 (see inset). Gray shaded areas depict the maximum 
possible range of MPA. (d) Typical time-of-flight images of the 
driven ID lattice for positive tunneling (top panel), negative 
tunneling (middle panel), and a heated, incoherent sample 
(bottom panel). 


£i{q) ~ El — 2Ji cos(ag), with ei being the band-center 
energy. 

III. PERIODIC DRIVING 

Periodic driving of the system is induced by a sinu¬ 
soidal frequency modulation Su{t) = z/q sm{Qt) of one of 
the two laser beams. This modulation gives rise to a peri¬ 
odic motion of the potential along the lattice axis. In the 
CO- moving frame, the atoms experience an inertial force 
F(t) = Fq cos(Ot) Bx with an amplitude of Fq = MViv^a. 
Apart from trap and interactions, the neutral particles 
are described by the Hamiltonian 

^ 2 

^ + V\!,t{x) + xFo cosiVlt). (1) 

The driving term in Eq. (1) breaks the translational sym¬ 
metry of the lattice which can be restored by the gauge 
transformation described in Appendix A [34, 35]. In the 
new reference frame of the lattice the resulting Hamilto¬ 
nian can be written as 

= ( 2 ) 


states |(a g), this yields a tight-binding Hamiltonian 

“ A{t)/h)\a q){a q\ 

“ V- 1 

FaFocos{ftt)^r]a,a'\c^'Q){c^ q\ , 

Oi' 

Here, and A{f) = and r]a,a' are dimensionless 

dipole matrix elements coupling bands a and a' defined 
in Appendix A. 

For driving frequencies that are large compared to the 
widths of the bands, the impact of periodic forcing can be 
understood as a combination of two effects as detailed in 
Appendix B. The first effect is a modification of each sin¬ 
gle band, described by the time-averaged single-particle 
dispersion relation 

^ “ ^( 0 /^) ( 4 ) 

= ea - 2JaJo{K) cos{aq). (5) 

This yields the effective modification of nearest-neighbor 
tunneling matrix elements = JasJo{K)^ where Jo 
is the zeroth order Bessel function of the first kind and 
K = aFo/(Ml) denotes a dimensionless driving ampli¬ 
tude. The second effect, resulting from the second term 
in Eq. (3), can be understood as the resonant coupling of 
the so-modified bands. 

Experimentally, the driving amplitude K is linearly 
increased to its final value within 50 ms after the atomic 
ensemble is prepared in the lowest-energy band of the op¬ 
tical lattice. Driving is maintained at the final amplitude 
for another 20 ms. Subsequently, all trapping potentials 
are rapidly switched off, and atoms fall freely under the 
influence of gravity for 40 ms time-of-flight before a res¬ 
onant absorption image is taken. 

For the lattice depths used throughout the presented 
experiments the atomic ensemble remains in the weakly 
interacting superfluid regime. Thus, the inversion of the 
effective band structure due to the sign change of for 
sufflciently large forcing amplitudes [see Fig. 1(b)] can be 
identified by the position of the coherent quasimomentum 
peaks in time-of-flight absorption images as shown in the 
upper two panels in Fig. 1(d). 

IV. INTERBAND MULTIPHOTON 
TRANSITIONS 

As explained in the previous section, beyond the tun¬ 
neling renormalization, the periodic forcing of frequency 
Q also induces finite matrix elements for coherent inter¬ 
band coupling processes that conserve quasimomentum 
but allow for energy changes of integer multiples of the 
“photon” energy hft [36]. Hence, an n^^-order multipho¬ 
ton transition is expected to occur when the resonance 
condition 


where the effect of the periodic driving is incorporated 
into a time-dependent vector potential A(t) given via the 
relation F(t) = —dA{t)/dt. In the basis of static Bloch 


1 AEf{q = ±7r/a) for jf < 0 


(6) 
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(a) }f=+OAia Ja^=0 Ja^=-0.4Ja /'^=+0.3 Jq 



Lattice depth Vq (units of Erec) 


FIG. 2. (Color online) Systematic investigation of multipho¬ 
ton spectra in the driven ID lattice, (a) The effective tunnel¬ 
ing matrix element acquires different values at the four 
measured driving amplitudes K, scaling with the Bessel func¬ 
tion J'o(K). (b) Excitation spectra for the four driving am¬ 
plitudes depicted in (a) with the maximum optical density 
encoded in brightness. Solid black lines, numbered on the 
right hand side, indicate the calculated positions of MPA to 
the first excited energy band according to Eq. (6) . 


with AE^{q) = e^{q) — £o^{q) is fulfilled. Condition 
(6) is plotted in Fig. 1(c) for a = 1 and n = 1,..., 10 
(white lines) as a function of the lattice depth for a driv¬ 
ing amplitude of K = 3.82 where = —0.4 Jq, (see 
inset). 

These interband transitions, induced by the periodic 
driving, significantly reduce the maximum optical den¬ 
sity extracted from the time-of-fiight images as can be 
observed in the bottom panel of Fig. 1(d) [36]. Such a 
decrease in the visibility in the absorption images can be 
ascribed to two distinct processes. First, further MPA 
might populate higher lying bands that are no longer 
trapped in the optical lattice and, thus, result in a de¬ 
crease of the optical density due to atomic losses. Second, 
interacting Bose-Einstein condensates in excited bands 
rapidly decay due to scattering processes, thereby reduc¬ 
ing the degree of coherence in the system [37, 38]. 

A spectroscopic study of these multiphoton transitions 
is shown in Fig. 2 where excitation spectra are obtained 
for various lattice depths at four different driving ampli¬ 
tudes. The resulting effective tunneling parameters 
[see Fig. 2(a)] yield different transition energies accord¬ 
ing to Eq. (6) that are plotted as solid black lines to¬ 
gether with the excitation spectra in Fig. 2(b). These ab 
initio calculations of transition energies exhibit excellent 
agreement with the experimental data. With increas¬ 
ing driving amplitude, higher-order excitations appear in 
the spectrum and the width of the resonances increases. 
While no excitations are present in the system above the 
third-order resonance for K = 1.7, up to the ninth or¬ 
der multiphoton transition can be identified at a driving 
amplitude of K = 7.0, where the lower orders of MPA 



(b) Driving frequency 0/2n (kHz) 



Driving frequency 0/2n (kHz) 



Driving frequency 0/2n (kHz) 


FIG. 3. (Golor online) Emergence of multiphoton interband 
transitions for increasing driving amplitude, (a) Multiphoton 
excitation spectrum obtained at a fixed ID lattice depth of 
9.5 Erec for increasing driving amplitudes. The positions and 
maximum possible widths of the depicted n^^-order transi¬ 
tions are indicated above the spectrum. Data a,t K ^ 4.0 
(second row from top) have not been measured and are in¬ 
terpolated. (b) Numerical simulation of the observed multi¬ 
photon interband excitation spectra. Plotted is the minimum 
occupation Nq of the lowest band observed during 20 ms of 
driving at the final value of K. Avoided crossings and exci¬ 
tations to higher energy bands are clearly visible. In (c) the 
resolution of the simulation data shown in (b) has been re¬ 
duced to match the resolution of the experimental data. To 
take into account the linear ramping to the final driving am¬ 
plitude K, excitations present for smaller values of K are kept 
in the spectra for larger K. The resulting spectrum clearly 
matches the experimental data depicted in (a). 


transitions already overlap significantly. One can also ob¬ 
serve that resonance features become weaker for increas¬ 
ing lattice depth; this is a consequence of the fact that 
the coupling matrix elements r]a,a' become smaller for 
deeper lattices. Each of the four individual data sets has 
been normalized, yielding comparable results for all in¬ 
vestigated driving amplitudes. Even for K = 2.4, where 
the effective tunneling amplitude vanishes, the signal-to- 
noise ratio is sufficient to clearly identify resonance fea¬ 
tures despite a significantly reduced level of coherence. 

In order to gain a deeper understanding of the emer¬ 
gence of interband MPA processes, the observed exci¬ 
tations are further explored with respect to the driving 
amplitude in Fig. 3. The excitation spectra depicted in 
Fig. 3(a) are obtained at a fixed one-dimensional (ID) lat¬ 
tice depth of Vb = 9.5 Erec while the driving amplitude K 
is gradually increased. Here, we observe that for increas- 
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ing final coupling strengths, the widths of the observed 
resonance features increase. Small additional features of 
the excitation spectrum in the region of large driving am¬ 
plitudes K and frequency Vt that cannot be explained by 
the expected resonance positions indicate the occurrence 
of multiphoton excitations to even higher-energy bands. 

A numerical simulation of the interband excitation 
spectra, described in more detail in Appendix D, is de¬ 
picted in Fig. 3(b). Starting from an ensemble of states 
with quasimomenta distributed sharply around the min¬ 
imum of the effective dispersion relation at the respec¬ 
tive maximum value of the time evolution under the 
time dependent Hamiltonian given in Eq.(2) has been 
integrated over 20 ms. Excitations to higher bands are 
quantified by the minimum fraction of occupation A'o of 
the lowest band during the time evolution. The obtained 
spectra exhibit good agreement with the experimental 
data in Eig. 3(a), reproducing the observed multiphoton 
transitions to the first excited band with n = 3 to 7 
as well as additional resonance features that are associ¬ 
ated with multiphoton excitations to higher lying bands 
at larger driving frequencies. Here, the finite width of 
the momentum distribution (resulting from thermal fluc¬ 
tuations, interaction-induced quantum fluctuations, and 
the trap potential) plays a central role, since it is re¬ 
quired to explain the observed resonances of even pho¬ 
ton number between the two lowest bands (see Appendix 
C). Avoided-crossing-type structures, indicating the res¬ 
onant hybridization of the first excited band with higher- 
lying ones, can be identified in the simulated resonance 
features at the n = 3, 4, and 5 transitions. 

To facilitate the comparison of the rich spectrum ob¬ 
tained from numerical simulations and the experimen¬ 
tal data, we have reduced the resolution of the numeri¬ 
cal data to the experimental resolution in Eig. 3(c). To 
emulate the linear ramping of the driving amplitude to 
its respective final value, each row n of this plot in¬ 
cludes the combined product of all rows corresponding 
to smaller values of K according to min[A'o(A"n)]^ = 
]Xi=i With this, an excitation present dur¬ 

ing any point of the ramp remains in the system at any 
larger value of the driving amplitude K. The resulting 
spectrum clearly reproduces the features of the experi¬ 
mental data shown in Eig. 3(a) such that the observed 
broadening of the multiphoton resonances can be as¬ 
cribed to the linear ramping procedure of the driving 
amplitude K. 

V. INTERBAND EXCITATIONS IN THE 
TRIANGULAR LATTICE 

Eor the investigation of interband MPA processes in 
more complex lattice structures we extend our studies to 
a driven two-dimensional triangular lattice. As depicted 
in Eig. 4(a) the lattice is composed of three running wave 
laser beams of equal intensity intersecting in the x^-plane 
with linear out-of-plane polarizations [36]. Here, inertial 
forcing is induced by a sinusoidal frequency modulation 



Vy(kHz) Lattice depth Vq (units of E^ec) 


FIG. 4. (Color online) Driving in the triangular lattice, (a) 
Illustration of the experimental setup, lattice structure and 
elliptical forcing, (b) Tunneling renormalization of the diago¬ 
nal lattice bonds, which depend on the horizontal and vertical 
frequency modulation amplitudes h>x/y' The isotropic renor¬ 
malization condition of i>y — VSux is plotted as a dashed 
line [36] . (c) Excitation spectrum for a fixed isotropic driving 
amplitude of K = 3.82, indicated by the solid dot in (b). Ex¬ 
pected positions and widths for multiphoton transitions are 
plotted as solid black lines similar to Eig. 2. MPA resonances 
of up to the fifth order can be clearly identified. 

of two of the three laser beams, resulting in a periodic 
elliptical forcing of the rigid lattice potential. This al¬ 
lows adjusting the amplitude and sign of two effective 
tunneling directions, denoted and in Fig. 4(a), 
independently. The tunneling renormalization of J®^ is 
shown in Fig. 4(b) to depend on the horizontal and verti¬ 
cal frequency modulation components Ux and which 
determine the ratio of the major and minor axes of the 
elliptical forcing F(t) [36]. For the investigation of MPA 
in the driven triangular lattice, we focus on an isotropic 
forcing parameter of K = 3.82 along all lattice bonds. 
This corresponds to a negative effective tunneling of max¬ 
imal amplitude (J®^ = = —0.4 Jq,), which is central 

for the study of frustrated magnetism [13]. 

MPA resonances of up to the fourth order can be 
clearly identified in the excitation spectrum shown in 
Fig. 4(c) with excellent agreement with the calculated 
transition energies to the first excited Bloch band (solid 
black lines). In addition, faint remnants of a fifth order 
transition are visible below lattice depths of 4.5 Erec- The 
maximum coherence of the atomic ensemble is reached 
only at a narrow parameter region between the fourth 
and the fifth order transitions while the maximum opti¬ 
cal density remains small for smaller driving frequencies 
where heating due to the resonant creation of collective 
intraband excitations might occur [29] . Such a significant 
limitation of the accessible parameter space for coherent 
manipulation of atomic ensembles is a crucial aspect for 
the experimental realization of periodic driving schemes 
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in two- or higher-dimensional lattice systems that rely 
on the applicability of time-averaged effective models. In 
addition, avoiding possible MPA processes is even more 
demanding for driving schemes employing more than a 
single driving frequency [39]. 


VI. CONCLUSION 

To conclude, multiphoton interband excitations have 
been investigated systematically with ultracold quantum 
gases in optical lattices. Thereby, multiphoton tran¬ 
sitions to the first excited energy band of up to the 
ninth order could be observed in a driven one-dimensional 
lattice as well as in a two-dimensional triangular lat¬ 
tice. The resonance positions are found to be in excel¬ 
lent agreement with ab initio calculations of the time- 
averaged effective single-particle band structure. Also, 
the strength of the resonances and their dependence on 
the driving amplitude show good agreement with numer¬ 
ical simulations. 

Our findings provide essential insights concerning the 
applicability of strong driving schemes for the experimen¬ 
tal realization of exotic quantum phases in the rapidly 
growing field of Floquet engineering. Moreover, a com¬ 
prehensive understanding of driven mesoscopic matter 
waves is a crucial prerequisite for the coherent control 
and addressability of intriguing quantum states. For in¬ 
stance, in analogy to coherent light-matter interactions, 
external periodic driving with precisely adjusted pulse 
shapes could allow for the generation of cat-like states 
between different Bloch bands in optical lattices [40] . In 
addition, the particularly strong inertial forcing enables 
the emulation of extremely strong field conditions in con¬ 
densed matter systems that are hardly accessible with 
real solids [33]. A further strength of quantum gases re¬ 
lies on the precise control over the interactions in the 
system. It is essential for the investigation of the com¬ 
plex interplay between periodic driving and interactions, 
which has very recently triggered several theoretical stud¬ 
ies [41-45]. 
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Appendix A: The driven lattice system 

In the following we provide a theory of multiphoton 
interband transitions in the driven ID lattice and will 
identify the basic processes that lead to multi-“photon” 
interband transitions in a driven optical lattice and their 
rates. 


In the lattice frame of reference a sinusoidal shaking 
of the lattice gives rise to a time-periodic homogeneous 
inertial force, resembling an ac voltage. Thus, a particle 
of mass M in a shaking cosine lattice is described by 
Eq. (1) of the main text: 

Ti{t) = ~ Y cos(27ra:/a) + xFq cos{Q,t), (Al) 


with Vb being the lattice depth, a being the lattice con¬ 
stant, and Fq being the amplitude of the force. 

It is convenient to describe the driven lattice in terms 
of the maximally localized Wannier states \ia) of the un¬ 
driven Hamiltonian {Fq = 0). The integer £ labels the 
lattice minima xi = £a and the index a = 0,1, 2,... de¬ 
notes the Blochs bands spanned by the corresponding 
Wannier states, with the energy increasing with a. The 
Wannier wave functions Wa{x — £a) = {x\£a) are real, ex¬ 
ponentially localized on a length increasing with a, sym¬ 
metric (antisymmetric) for even (odd) bands, Wa{—x) = 
{—)^Wa{x)^ and shall be normalized, J dx \ wa{x)\‘^ = 1. 

In Wannier representation, non-interacting particles 
in the driven lattice are described by the tight-binding 
Hamiltonian 


m = E ea|a£)(a^| - Y (|q;(^ + 2))(q;^| + h.c.) 

(^\at){a£\ + £){a£^ 


£a 

- K cos( 


(A2) 


The Hamiltonian is characterized by the amplitude of the 
potential modulations K = Foa, the band-center energies 


/ 


dx Wa{x) 


2M 


,di 


Vq 


—^cos(27rx/a) Wa{x)^ 


the nearest-neighbor tunneling parameters 


-/ 


Ja = — dx Wa {x — a) 


2M 


di 


Vq 


—^cos(27rx/a) Wa{x), 


and the dimensionless dipole matrix elements 


(A3) 


(A4) 


Va'a = - [ dxWa'{x)xWa{x) (A5) 

a J 


describing interband coupling. As a consequence of the 
parity of the Wannier functions, r]a'a is non zero only 
if {a' — a) is odd. Moreover the sign of the tunneling 
parameter alternates with the band index Ja/\Ja\ = 
(—)^. It is an approximation to neglect tunneling beyond 
nearest-neighbor sites and not to take into account time- 
periodic band-coupling terms between Wannier states on 
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different lattice sites. However, for the two low-lying 
bands those terms are rather small. The neglected terms 
might still be relevant for the description of higher-lying 
bands though. Including the neglected terms in the anal¬ 
ysis presented below would be straightforward, but is not 
done here for the sake of a simple presentation, captur¬ 
ing the basic picture. However, these terms are included 
in our numerical simulation of interband transitions de¬ 
scribed in Appendix D. 

The driving term breaks the translational symmetry of 
the lattice. However, the symmetry can be restored by 
performing a gauge transformation 


u’{t) = U\t)H{t)U{t) - ihU\t)U{t), 
with the unitary operator given by 

U (t) = exp 

£a 


and 


K£ 

X^(t) = —— / dt'cos{ftt') =—iK 

d Jo 


(A6) 

(A7) 


(A8) 


where K = K/(Ml). This unitary transformation inte¬ 
grates out the oscillatory shift in quasimomentum by 

Ag(t) = — — sin(f^t) (A9) 

a 

induced by the periodic force. The new Hamiltonian 
reads 


n'it) = 


Sa\o:i){ai\ — + l)){a£\ + H.c.) 


ia *- 

a' 


(AlO) 


with time-periodic Peierls phase 0{t) = x^(t) = 

—aAq{t) = K sin(Ot). 

Let us, for simplicity, assume a system of N lattice 
sites under periodic boundary conditions and express the 
Hamiltonian in terms of Bloch states \aq), with the quasi¬ 
momentum quantum number q taking N discrete values 
k = 27r/i/(Aa), with integer /i in the interval ( —y, y]- 
Using {a'£\aq) = exp(t^ag), one finds 


= (All) 

q 


with 




ea {q - Mt)/fi)\<^q){<^q\ 


+ Kcos{Q.t) y] Va'aWq){aq\ 


(A12) 


corresponding to Eq. (3) in the main text. Here, 

A{t) = —hAq{t) (A13) 


plays the role of a vector potential and 

= £a — ‘^Ja cos(ga) (^14) 

denotes the single-particle dispersion relation for band a. 

As a consequence of the discrete translational invari¬ 
ance, the Hamiltonian {t) conserves quasimomentum 
q. That means that interband transitions will occur be¬ 
tween Bloch states \aq) and \a'q) of the same quasimo¬ 
mentum q. This reduces the problem to independent 
subspaces that are characterized by q and spanned by 
states labeled by the band index a. 


Appendix B: Floquet picture 

Let us investigate the problem within the extended 
Floquet Hilbert space. The time-dependent Schrodinger 
equation 

iMtm)) = 'fL'{t)m)) (Bi) 

possesses a set of generalized stationary states of the form 

IV’r.(i)) = (B2) 

called Floquet states. They are characterized by the 
time-periodic Floquet mode \uy{t)) = \uy{t -\-T)) and by 
the quasienergy 5^^, where the driving period T is defined 
by T = 21 :jVt. The Floquet states, which are labeled by 
some quantum number form a complete basis of the 
state space at any time t. Therefore, we can expand the 
time evolution of a state in terms of the Floquet 

states like 

IV’(i)) = E = E (B3) 

U V 

with time-independent coefficients = ('0zy(to)|'0(to)). 

Unlike the Floquet states, the quasienergies and the 
Floquet modes are not defined uniquely. For each 
Floquet state a whole family of Floquet modes and 
quasienergies, 

evm = £!/ + mWt, \uum{t)) = (B4) 

labeled by the integer m, can be defined such that 

(B5) 

for all m. However, Floquet modes and quasienergies of 
different m still constitute independent solutions of the 
eigenvalue problem 

Q\'^vrrif) — ^vrn\'^vrrifl 
of the quasienergy operator 

(B7) 

This eigenvalue problem is defined in the extended Flo¬ 
quet Hilbert space, being the product space of the state 
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space with the space of time-periodic functions. In this 
space the scalar product is given by 

{{u\v)) = ^ dt{u{t)\v{t)). (B8) 

When considering a periodically time dependent state 
\u{t)) = \u{t^T)) as an element of the extended Hilbert 
space, we write it as a double-ket Likewise, an op¬ 
erator is marked by an overbar, like Q, if it is considered 
to act in the extended Hilbert space. 

For the driven lattice a useful set of basis states span¬ 
ning the extended Hilbert space is given by 

\aqm)) : \aqm{t)) = \aq)e^'^^^ (B9) 

and labeled by the band index a, the quasimomentum q, 
and the Fourier index m. These states are the Floquet 
modes of the undriven problem with K = 0. With re¬ 
spect to these basis states the quasienergy operator pos¬ 
sesses the matrix elements 


{{a'q'm'\Q\aqm)) = {a'q'\{'H'^,_^ +6m'mfi^)\aq), 

(BIO) 

where 


n'^ = ^ (Bii) 

denotes the Fourier transform of the Hamiltonian such 
that H'{t) = The quasienergy operator 

assumes a block structure with respect to the index m, 
which plays the role of a relative photon number. The 
diagonal blocks describe subspaces of different “photon” 
number and are shifted relative to each other in energy by 
integer multiples of the photon energy h£L. The diagonal 
blocks are coupled by off-diagonal blocks characterized 
by describing m-photon processes. 

Let us evaluate the matrix elements explicitly. One 
finds 


= ^^<xm{q)\aq){aq\ + \K5\rn\,i'ncc'<x\a'q){aq\ 


qa 


where 


(B12) 


1 y 

Samiq) = 7f dtea{q- A{t)/h), 

^ Jo 

= ScSmfi - JaJmiK) 

(B13) 

with Jjn denoting a Bessel function of order m. The 
diagonal blocks are given by 

{{a'q'm\Q\aqm)) = 5q>qSa>a [e®®(g) + mhO], (B14) 


where we have introduced the effective dispersion relation 
^‘fiq) = £ao{q) = £« - 2JaJo{K) cos{aq). (B15) 


IlgO)) 


\0q2)) 


nhQ. 
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^£eff(<?) 
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\'i-q-n}) 

\lqi-l-n))) 


ha 


FIG. 5. (Color online) Resonance condition for an n-photon 
transition between the bands o = 0 and o = 1 in the extended 
Floquet Hilbert space. The energy levels correspond to the 
unperturbed quasienergies e^{q) -h mhfl. The states \aqm)) 
are labeled by the band index a, the quasimomentum wave 
number q, and the relative “photon” number m. 


The off-diagonal blocks read 


{{a'q'm'\Q\aqm)) 


= S, 


^a'a^a{m' —m) {o) 


T 2^^^\m' —m\,lVo!.'( 


(B16) 


The first term describes multi-photon processes with ar¬ 
bitrary \m' — m\ that do not change the band index a. 
The second term describes single-photon transitions with 
|m' — m| = 1 between two bands a and a' with odd 
{a' — (a). Thus, the quasienergy operator does not di¬ 
rectly possess matrix elements that describe multiphoton 
interband transitions. Such processes emerge, however, 
from-higher order virtual processes in the extended Flo¬ 
quet Hilbert space, as we will discuss in Appendix C. 


Appendix C: Multiphoton interband coupling 

The basis states \aqm)) correspond to eigenstates \aq) 
of the undriven Hamiltonian. An m-photon interband 
coupling process from state \aq) to state \a'k) is expected 
to occur when the resonance condition 

- ef{q) =nhn + 5 (Cl) 

with integer n and sufficiently small detuning 5 is ful¬ 
filled. In the extended Floquet Hilbert space this reso¬ 
nance condition corresponds to a quasidegeneracy of the 
unperturbed states \aqm)) and \a'q{m — n))) with respect 
to the diagonal blocks. The energy cost of the tran¬ 
sition is compensated by the destruction of n photons. 
This is illustrated in Fig. 5 for the case of a = 0 and 
a' = 1, which is relevant for the experiment. Resonant 
n-photon interband excitations are expected when the 
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FIG. 6. (Color online) Single-photon transition. 


(effective) coupling matrix element between both states 
becomes comparable to the detuning 6. If both states 
are not coupled directly by a matrix element appear¬ 
ing in the off-diagonal blocks, they might be coupled via 
an effective matrix element resulting from higher-order 
processes via energetically distant “virtual” intermediate 
states \a''qm'')). 

Let us estimate the relevant coupling matrix element 
for an n-photon process from the lowest into the first ex¬ 
cited band at quasimomentum q, i.e., between the states 
|0g) and |lg). For this purpose, we have to evaluate 
the coupling matrix element between states | 0 g 0 )) and 
[Iq — n)) in the extended Floquet Hilbert space. For a 
single-photon process with n = 1 we find a direct cou¬ 
pling matrix element 


C^ = {{lq-l\Q\OqO)) = JfK (C 2 ) 


as illustrated in Fig. 6 . 

For a two-photon process with n = 2, we do not find a 
direct coupling matrix element. However, both states can 
be coupled via two second-order processes, \0q0)) —> | 0 g — 
1 )) ^ |lg — 2 )) and |0g0)) ^ |lg — 1 )) ^ |lg — 2 )) as shown 
in Fig. 7. Using the standard expression of degenerate 
perturbation theory the effective coupling matrix element 
is given by 


C2 = {{lq-2\Q^^\0q0)) 
£o-i{q) 


= 


Ml 


nn )' 


(C3) 


neglecting 5 in the energy denominator (this level of ap¬ 
proximation is equivalent to a high-frequency approxima¬ 
tion [29]). The order of magnitude of this term can be 
estimated by noting that for small arguments the Bessel 
function behaves like Jm{x) ^ such that for m 7 ^ 0 


/ ^ \ 1^1 


Including also the momentum dependence of 
which for odd m is just given by a factor of sm{aq), we 
find that the matrix element of the two-photon process 


FIG. 7. (Color online) Two-photon transition. 



FIG. 8. (Color online) Three-photon transition. 


is of the order of 

C 2 ^ sm{aq) Ji, 2 . (C5) 

Typical paths contributing to three-photon interband 
transitions in leading order are depicted in Fig. 8 . They 
also involve intermediate states of higher-lying bands and 
give rise to effective tunneling matrix elements 

C 3 = ((19 - | 0 g 0 )) ~ K. (C 6 ) 


Generally, the coupling matrix elements describing an 
n-photon interband transition obey 

n —1 

k for odd n (C7) 




nn 


and 


Cn ^ sin(ag) ( ^ for even n. (C 8 ) 
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Driving frequency Q/2n (kHz) 


FIG. 9. (Color online) Minimum occupation of the lowest 
band during 20 ms of time evolution, versus driving frequency 
Q, and driving strength K. At t = 0 the system is prepared in 
the Bloch state corresponding to the minimum of the effective 
dispersion relation of the lowest band, that is with quasimo¬ 
mentum q = 0 (g = Tv/a) for K lower (larger) than 2.4. 


(a) 



(b) Driving frequency Q/2n (kHz) 



Driving frequency Q/2n (kHz) 


The factor of sin(ag) results from the fact that for even 
n always one matrix element Sam{Q) with odd m always 
contributes in leading order. [Also higher-order coupling 
paths for even n contain at least one factor SamiQ) oc 
sin(ag) with odd m.] 

We can see that transitions involving an even num¬ 
ber of photons are suppressed by an additional factor of 
sm{aq)Ja/{hQ) with respect to transitions with odd pho¬ 
ton numbers. In particular for the experimentally rele¬ 
vant quasimomenta q = 0 and q = nja resonances with 
even n are suppressed completely. This is a consequence 
of the fact that only bands with Wannier functions of 
different parity are coupled directly by the periodic force 
on the level of our approximation. The transitions with 
even photon numbers n are observed experimentally be¬ 
cause of the broadening of the quasimomentum distri¬ 
bution due to temperature, interactions, and the finite 
system extent. 


Appendix D: Simulation of the time evolution 


In order to integrate the time evolution of the shaken 
lattice, let us start directly from the Hamiltonian in real- 
space representation (Al). As before, we perform a gauge 
transformation to restore the translational symmetry of 
the lattice 

n'{t) = u{tyn{t)u{t) - ihu\t)u{t) (di) 


with 


U{t) = exp (- - 


(— t y dt'F qx cos{nt)^ 


= exp 


(^iAq{t)xJ 

(D2) 


giving 

where A{t) again plays the role of the vector potential. 
Assuming a system of length L with periodic boundary 


FIG. 10. (Golor online) Same as Fig. 9, but with the initial 
state shifted away from the minimum by (a) Aq = O.Itt and 
(b) Aq = 0.27r. 


conditions, we can express the Hamiltonian in terms of 
momentum eigenstates \p) with wave functions 

x\p) = ^exp{ipx). (D4) 

For that purpose it is convenient to decompose the mo¬ 
mentum wave number as 

p = q + PQ (D5) 

with — ^ < q < ^5 being an integer, and Q = 27r/a. 
Introducing the localization energy F^ioc = j (2Ma^), 

which describes the kinetic-energy cost of localizing a 
particle on a lattice constant a, as the natural unit of 
energy, we find matrix elements 

(g' + /?'Q|7i'(i)|g + /3Q) = 5qi^qhi3ip{q,t)Eioc, (D6) 

with 

hi}'0{q,t) = 50:i}^[q +PQ - A{t)/h\^ 

+ ^(V,/3+i + (D7) 

One can see that the wave number q is conserved, so 
that the dynamics occurs in the space spanned by the 
integer quantum number By diagonalizing the dimen¬ 
sionless Hamiltonian h^'^i3{qA) for A{t) = 0, we obtain 
the band structure of the undriven lattice. The fact that 
the diagonal matrix elements increase like 4/3^, while 
the off-diagonal terms are constant, = Vo/(4^ioc)5 

shows that Bloch states with energy much larger than the 
lattice depth resemble plane waves. 

For our simulation we initialize the system in a Bloch 
state of the lowest band with quasimomentum q and in¬ 
tegrate the time-dependent Schrodinger equation of the 
dimensionless time-dependent Hamiltonian hp'p{qA) for 
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0 2 4 6 8 10 12 14 16 18 20 


Driving tinne (nns) 

FIG. 11. (Color online) Occupation of the lowest bands during 
the time evolution. Plotted is the overlap squared of the time 
evolved state ^{t) with eigenstates (j)j for (a) K — 3.0 and 
Q = 27rx4.65kHz, (b) K = 3.0 and Q, = 27rx5.00kHz, and 
(c) K — 2.5 and = 27rx2.90 kHz. The bands for a = 0 to 
3 are depicted with increasing brightness (see legend). Bands 
above a = 3 exhibit no signihcant occupation and are omitted 
from the plots. 


the given forcing strength K over a time span correspond¬ 
ing to 20 ms. We take into account 61 plane waves; in¬ 
creasing this number further does not alter the results. 
We report the minimum occupation of the lowest band 
encountered during the time evolution. This number pro¬ 
vides a measure for the amount of interband excitation 
expected on the given time scale. 

We choose q to be the minimum of the effective disper¬ 
sion relation of the lowest band, that is q = 0 {q = tt/ a) 
for K less (greater) than 2.4. With respect to the driving 
frequency and strength, the minimum occupation of the 
lowest band is plotted in Fig. 9. While this plot already 
resembles the experimentally measured data in some re¬ 
spects, it hardly shows resonances corresponding to even 
photon numbers n. This suppression of even resonances 
is expected from the theory presented in Appendix C. 

The even resonances observed in the experiment can 
be explained by the finite width of the quasimomentum 


distribution, induced by the finite extent of the trapped 
system, finite temperature, and interactions. In order to 
take into account that the quasimomentum distribution 
possesses a finite width re, we simulate the time evolution 
also starting from initial states that are shifted away from 
the minimum of the effective dispersion relation by Aq. 
As examples, data for Aq = 0.l7r/a and Aq = 
are shown in Figs. 10(a) and 10(b), respectively. Here 
even resonances are clearly visible. The plot shown in 
Fig. 3(b) is a superposition of resonance data obtained 
for different values of Aq (varied in steps of 0.0257r/a), 
with Gaussian weights oc exp(—Ag^/re^). 

The width of the momentum distribution was set to 
w = 0.l7r/a. While this value cannot be determined 
experimentally with sufficient accuracy since the time- 
of-flight pictures shown in Fig. 1(d) are taken before the 
far-held limit is reached [46], Fig. 1(d) still provides an 
upper bound for re, which is of the order of 0.l7r/a. An 
estimate for the lower bound of w is obtained from the 
Thomas-Fermi radius and is of the order of 0.027r/a [36]. 
Since thermal and quantum huctuations will cause fur¬ 
ther broadening, the value w = 0.l7r/a is a reasonable 
assumption. 

In Fig. 9 we can see that the three-photon resonance 
(near 5 kHz) from the ground band a = 0 into the hrst 
excited band {a = 1) is split into two resonances. This is 
a signature of the fact that the hrst excited band is cou¬ 
pled resonantly to even higher lying bands, so that an 
avoided crossing is formed in the quasienergy spectrum. 
The doubled resonance rehects this avoided crossing and 
explains the large broadening of the three-photon reso¬ 
nance visible in the experimental data [Fig. 3(a), main 
text]. In order to identify the bands involved in the 
three-photon transition, let us have a look at the sim¬ 
ulated time evolution. In Fig. 11 we plot the occupations 
of the lowest bands over time for diherent parameters K 
and ft for an initial state with q = nja. At the driving 
strength K = 3 [see Fig. ll(a,b)], the left- and right-hand 
sides of the double resonance are captured roughly by 
Cl = 27rx4.65kHz and Cl = 27rx5.0kHz, respectively. We 
can see that the system performs Rabi-type oscillations 
between the ground band and a hybridized state with 
strong contributions from the three bands with a = 1, 2 
and 3. We can compare these results with a plot of the 
time evolution near the five-photon resonance {K = 2.5, 
Cl = 271 x2.9kHz) shown in Fig. 11(c), where only the 
first excited band becomes populated significantly, since 
for this resonance an avoided crossings occurs only for 
larger K (Fig. 9). 
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Supplemental Material 


The following supplemental material provides additional information on the experimental 
setup and preparation procedure as well as the employed time-periodic driving scheme 
and data evaluation. 


SI. EXPERIMENTAL SETUP 

All experiments presented in this paper start with a 
Bose-Einstein condensate (BEC) of atoms in an el¬ 
liptical crossed optical dipole trap operating at a wave¬ 
length of Adt = 1064 nm. With minimal beam waists 
of wo,h = 245 fim horizontally and wo,v = 82 fim verti¬ 
cally, final trapping frequencies of cjh = 27rxl9Hz and 
ujy = 27rx48Hz are reached after evaporative cooling to 
quantum degeneracy. Eollowing the creation of a BEC 
consisting of (3±1) x 10^ atoms in the crossed dipole trap, 
the ID optical lattice intensity is linearly ramped up to 
its final value in 100 ms. With a wavelength of at = 
830 nm the recoil energy of the lattice amounts to lErec = 
{2M) = hx 3.33 kHz, where M = 1.44x10“^^ kg 
denotes the atomic mass of ^^Rb. As the two running- 
wave laser beams of the ID lattice enclose an angle of 
'd = (117.1±0.2)° the resulting lattice constant is given 
by a = Al /[2 sin('d/ 2 )] = (486.5 ± 0.5) nm. Eor ^^Rb 
atoms in the |E= 1 , mF= — l) hyperfine ground state with 
an inter-particle scattering length of = (100.4 ±0.1)ao 
the Thomas-Eermi profile [47] for the given particle num¬ 
bers and trapping frequencies yields a number of 59- 
69 occupied lattice sites in pancake-like shapes with in¬ 
creasing particle numbers at the center of the trap (see 
Tab. SI). Hereby, the occupation-weighted on-site inter¬ 
action Uw = UiTii ranges from 1.33xlO“^Erec to 

1.01 X10 ^E^rec- Similar experimental conditions are used 
in Ref. [15]. 


S2. TRIANGULAR LATTICE 

The triangular lattice [13, 48] is comprised of the two 
laser beams that are also used for the creation of the 
ID lattice (denoted beam 1 and 2 in the following) with 
an additional running-wave laser beam (denoted beam 


3) along the y axis. Ideally, all three beams in the xy- 
plane should be aligned at 120 ° with respect to each 
other, yielding wave vectors of ki /2 = ^l/2 (t\/3, “I 5 O) 
and ks = /cl (0,1,0). With each beam being linearly 
polarized perpendicular to the lattice plane, this setup 
results in a triangular lattice with a lattice spacing of 
% = 2 Al/ 3 = 553.3 nm. 

The resulting scalar lattice potential reads 


V(r) 


Vo 


3 

6 + 4y]cos(bjr-A(/)jfc) , 

i=l 


(SI) 


where are the phase differences between the laser 

beams with wave vectors kj and k/^. The denote three 
reciprocal lattice vectors that are given by the circular 
permutation 


bz — ^ijk (kjf ^A:) ( 82 ) 

where £ijk is the Levi-Civita symbol, such that bi = 
h ( 1 , 0 , 0 ) and b 2/3 = 5/2 (—1, =F\/3, O) with h = VSkL^ 
Note that the Einstein summation convention is not ap¬ 
plied in Eq. (S2). Bravais lattice vectors in real space 
can simply be obtained from pairs of reciprocal lattice 
vectors with the relation • hj = 27t 6ij such that, e.g., 
ai = % (0, —1, 0) and a 2 = a^/2 (\/3, —1, O). 

However, as mentioned above, the alignment of the 
lattice beams is slightly distorted due to a limited opti¬ 
cal access in the experimental setup. The experimentally 
realized angles between the three laser beams can be de¬ 
termined from the position of quasimomentum peaks in 
time-of-flight images that are generated by a Kapitza- 
Dirac diffraction and are given in Tab. S2 together with 
the corresponding lattice constants. Eor both the ID lat¬ 
tice as well as the triangular lattice, all calculations of 
band gap energies shown in the main text are performed 
with the exact measured lattice beam angles. 


A 

2x10® 

3x10® 

4x10® 

Wave vectors 

Angles 'dij 

Lattice constants 

A Sites 

59 

65 

69 

/(ki,k2) 

117.1° ±0.2° 

a 2 = 545.9 ± 0.8 nm 

U/max 

6.3x10® 

8.7x10® 

1.1x10"^ 

x:(ki,k3) 

120.4° ±0.4° 

|ai| = 560.6 ± 1.0 nm 

Uw 

1.33xlO“^Srec 

l.lSxlO-^Srec 

1.01xl0“‘‘Srec 

/(k2,k3) 

122.4° ±0.4° 

ai — a 2 = 554.3 ± 1.0 nm 


TABLE SI. System parameters of the one-dimensional lat¬ 
tice with number of atoms A, number of occupied lattice 
sites Asites, maximum site occupation Umax and occupation 
weighted on-site interaction parameter Uw 


TABLE S2. Measured angles between the lattice vectors k^ 
and the corresponding three lattice constants of the slightly 
distorted triangular lattice setup [compare Fig. 4(a) in the 
main text]. 
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Preparation of the atomic ensemble in the triangular 
lattice follows the same procedure as for the ID lattice. 
Since the ensemble is only weakly confined along the 2 :- 
axis (perpendicular to the lattice plane) the atoms form 
an array of approximately 2100 — 2600 elongated cigar¬ 
shaped tubes. Similar experimental conditions are used 
in Ref. [17] . On average the tubes are occupied with 70-95 
atoms while the external harmonic confinement results in 
an occupation of 175-235 in the center of the trap which 
fades towards the edges of the ensemble. 


S3. SINUSOIDAL PERIODIC DRIVING 

In the ID optical lattice, periodic inertial forcing is 
induced by a frequency modulation v(t) = v-\-5v{t) with 
= z/q sin(Dt) of one of the two running-wave laser 
beams. The frequency modulation translates to a phase 
change of 


(j){t) = —27r f dt'6u{t') (SI) 

J—00 

and, thus, a global shift of the rigid lattice poten¬ 
tial according to Viat(0(t),r) = ldat( 05 r ~ I^(^)) with 
R(t) =ex(t){t)/b denoting the trajectory of a lattice well. 
Since F(t) = —MR(t) this leads to a periodic forcing of 

F(t) = —Fq cos (Dt) e^,, with Fq = MQuoa. (S2) 

This forcing scheme can be easily extended to the trian¬ 
gular lattice. Here, two of the three lattice laser beams 
(with wave vectors ki/ 2 ) are modulated according to 

^^ 1/2 (^) = sin (Dt) + I'y cos (Dt). (S3) 

The resulting elliptical forcing in the triangular lattice is 
given by 


F|>(t) = —Fx cos (Ut) Bx — Fy sin (Dt) By (S4) 

with independent forcing amplitudes along the x and y- 
axis Fx = V^MflUxCi and Fy = M^UxCi respectively. 
Hereby, the prefactor ^/3 in the expression for Fx stems 
from the different projection of the forcing onto the di¬ 
agonal lattice directions. Isotropic forcing along all three 
lattice bonds is, hence, reached if the frequency mod¬ 
ulation amplitudes satisfy the condition Uy = VSux as 
illustrated in Fig. 4(b) in the main text. 


S4. EXPERIMENTAL OBSERVATION OF 
TUNNELING RENORMALIZATION 
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FIG. SI. (Color online) Quasimomentum distribution of 
driven ID lattices, (a) Dispersion relation of a lattice with 
positive tunneling matrix element and corresponding mea¬ 
sured momentum distribution after time-of-flight for a super¬ 
fluid state, (b) Effective dispersion and measured momentum 
distribution for negative effective tunneling. The inversion of 
the dispersion relation is clearly visible in the absorption im¬ 
age by the absence of a zero-momentum component as atoms 
occupy the edges of the first Brillouin zone, (c) Definition 
of the momentum visibility mask for time-of-flight images of 
the atomic ensemble in the driven one-dimensional lattice for 
the characterization of the degree of coherence as well as the 
effective dispersion relation in the driven system. 


vectors results in a renormalization of the bare tunneling 
matrix elements Jij between adjacent lattice sites i and 
j according to 


Jtf = Jij ^ ^ di exp [iWij {t)/h] , (SI) 

where i denotes the imaginary unit and the projection of 
the forcing onto the lattice bonds is given by 

Wij{t) = - f dt'F{t')^Rij. (S2) 

J—00 

Since the tunneling renormalization is chosen to be real¬ 
valued and isotropic for all experiments presented within 
this manuscript, we write Jij = J and in 

the following. For sinusoidal forcing this yields the dis¬ 
cussed renormalization of = JJq {K). For a given 
lattice potential the dimensionless driving amplitude K 
depends only on the frequency modulation amplitude. 
In the case of the ID optical lattice with periodic forcing 
along the lattice axis according to Eq. (S2) all tunneling 
matrix elements renormalize with 


The renormalization of effective tunneling matrix ele¬ 
ments is essential for the interpretation of the obtained 
multiphoton spectra. As the presented considerations 
hold for all energy bands, we omit the band index a 
in the following. The periodic forcing onto the lattice 


aFo MaJ 


(S3) 


Hence, the frequency modulation applied here has the 
advantage of the tunneling renormalization Jo{K) be¬ 
ing independent of the chosen modulation frequency U 
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in contrast to inertial forcing schemes that rely on direct 
phase modulation as, e.g., realized by a piezo actuator 
attached to the retro-reflecting mirror of a counterprop- 
agating lattice setup. The renormalization of the tunnel¬ 
ing matrix elements with a zeroth-order Bessel function 
of the first kind implies that the effective tunneling be¬ 
comes negative for sufficiently large forcing parameters 
and, hence, the band structure is inverted. In Fig. SI 
both cases of positive- and negative effective tunneling 
matrix elements are presented [compare Fig. 1(d) in the 
main text]. Here, the time-of-flight (TOF) images clearly 
show a coherent occupation of the minima of the effective 
band structure. For positive tunneling [Fig. Sl(a)] the 
atomic ensemble exhibits a pronounced quasimomentum 
component at g' = 0 with integer multiples dX q = :^2'k ja 
in accordance with the corresponding Wannier envelope 
of the momentum distribution. In contrast, the inver¬ 
sion of the effective band structure results in the atomic 
ensemble residing at the edges of the first Brillouin zone 
dX q = ±.711 a for the case of negative effective tunneling 
without having a quasimomentum component of zero as 
depicted in Fig. Sl(b). It is important to note that the 
total driving time can be experimentally fine-adjusted 
such that the forcing function F(t) completes a full pe¬ 
riod in good approximation before all trapping potentials 
are shut off for the time-of-flight imaging. At this point 
in time, the momentum transfer of the forcing onto the 
atomic ensemble due to the micromotion vanishes such 
that the oscillating Wannier envelope coincides with the 
Wannier envelope of a system at rest [15]. 

A convenient quantification of the forcing-dependent 
behavior of the atomic ensemble is given by the visibility 
as defined in Fig. SI(c). The mask applied to the TOF 
images yields a positive (negative) contribution for an 
occupation of quasimomentum states corresponding to 
positive (negative) effective tunneling matrix elements. 
Hereby the summation scheme over the optical density 
image pixels ensures equal weighting of positive and neg¬ 
ative contributions. 

In Fig. S2 the tunneling renormalization is investigated 
with respect to the forcing strengths. Experiments are 
performed in a ID lattice with a lattice depth of Vb = 
S.OErec- Hereby, the bare tunneling amplitude in the 
lowest band amounts to Jq = 1.2x10“^ Erec- 

As the energy gap between the two lowest bands of the 
initial dispersion relation is Egap = 3.6 E^-ec = hx 12.1 kHz 
the driving frequency is chosen to Q = 27rxI.5kHz. 
TOF-images have been taken for an increasing final forc¬ 
ing parameter K. The column sum of these images is 
depicted in Fig. S2(a), showing a clear and sudden jump 
between the two distinct cases of positive and negative ef¬ 
fective tunneling. This jump is reproduced when extract¬ 
ing the momentum contrast from the images [Fig. S2(b)]. 
Zero crossings coincide extremely well with the ab initio 
calculation of the corresponding Bessel function Jq{K). 
Remarkably, the measurements still exhibit sharp super¬ 
fluid momentum peaks even after the third zero-crossing 
of the Bessel function, indicating that the coherence of 
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FIG. S2. (Color online) Periodic driving of the ID lattice 
with increasing amplitude, (a) Row-sum of measured TOF- 
images for increasing forcing parameter K. The jump in 
quasimomentum occupations at the sign-change of the ef¬ 
fective tunneling is clearly visible, (b) The zero-crossings of 
the extracted visibility (data points) coincide remarkably well 
with the zero-crossings of the renormalization Bessel function 
(dashed line) that are emphasized by dotted vertical lines. 
Error bars indicate the standard deviation of the obtained 
visibility data as each column is averaged over at least two 
individual TOF-images. 


the system is retained for very strong driving amplitudes 
as long as any multiphoton excitations can be avoided by 
a suitable choice of driving frequency and lattice depth. 
Exemplified, a forcing parameter of K = 10 corresponds 
to a driving amplitude in real-space of approximately five 
lattice sites. 

Eor the triangular lattice and an elliptical forcing given 
by Eq. (S4) the renormalization can be adjusted indepen¬ 
dently along the vertical lattice bonds and the two diag¬ 
onal bonds [denoted and compare Eig. 4(a)] by 
choosing the vertical and horizontal frequency modula¬ 
tion amplitudes Uy and Ux respectively. The correspond¬ 
ing forcing parameters are 

= -^Vy and + ^y- (§4) 

In Eig. 4(b) the renormalization for the diagonal tunnel¬ 
ing matrix elements is plotted in dependence of Ux 
and Uy. The isotropic inversion of all tunneling matrix 
elements results in a peculiar band structure with two 
degenerate minima at the vertices of the first Brillouin 
zone. Note that, due to the reduced total optical den¬ 
sity in the time-of-fiight images of the driven triangu¬ 
lar lattice, driving at the maximum amplitude following 
the 50 ms linear ramp is reduced to a duration of 2 ms. 
Properties of such a driven triangular lattice system are 
thoroughly discussed in the Refs. [13] and [17]. 
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FIG. S3. (Color online) Excitation spectra in dependence 
of the number Umax of maximum-valued pixels included for 
the calculation of the maximum optical density. Spectra for 
each value of included pixels are plotted as colored lines cor¬ 
responding to the color bar. The red solid lines indicate exci¬ 
tation spectra for a number of Umax = 10^ included pixels as 
used for Fig. 3(a) in the main text. 

S5. DATA EVALUATION 

While the visibility defined in the previous section is 
useful for the characterization of the effective tunneling, 
it cannot easily be applied to the investigation of multi- 
photon excitations in driven lattices. On the one hand 
this quantification method has two decisive advantages: 
To begin with, it yields per definitionem a normalized sig¬ 
nal due to the summation and normalization over equal 
surfaces of the absorption image. Furthermore, the ob¬ 
tained visibility value does not only reveal the level of co¬ 
herence in the system but also gives rise to a sign change 
of the signal if the effective dispersion relation is inverted 
in the time-averaged Floquet picture. On the other hand, 
it can be impractical to account for the periodic motion 
of the Wannier envelope in momentum space, e.g., by 
including more momentum components into the mask 
of Fig. S 1(c) or by fine-adjusting the driving times for 
varying driving frequencies. Moreover, the quality of the 
signal is prone to fluctuations. For example the location 
of the superfluid momentum peaks has to be precisely 
defined. Finally, the visibility does not yield meaning¬ 
ful results for the case of vanishing effective tunneling 
(compare Fig. 2 at A = 2.4 in the main text). 

For these reasons, the maximum optical density has 
been chosen to quantify the degree of coherence in the 


system. It is simply obtained by averaging the optical 
density values for a number nmax of maximum pixels in 
each TOF image: 

-| '^max 

Max.ODj =- W Pi (SI) 

^max . 

1=1 

with the set of Up ^ rimax optical density values for each 
pixel P/ = {^ 1 ,^ 2 , • • • :Pnp} that are sorted in descending 
order for each TOF image I. Although this method does 
not yield any information concerning the momentum of 
the atomic ensemble, it is, thus, an excellent tool for 
measuring excitation processes for systems that result in 
an overall loss of coherence. 

The number of included pixels that are averaged to ob¬ 
tain the maximum-optical density is arbitrary. A small 
value of Umax may result in a signal that is very sen¬ 
sitive to even small excitations in the investigated sys¬ 
tem but can also be prone to detrimental fluctuations 
such as faulty pixels of the CCD. In addition, the to¬ 
tal value of the maximum optical density can also be 
strongly influenced by the renormalization of tunneling 
in the effective time-averaged Floquet picture for small 
^max- While for the excitation spectra shown in Fig. 2(b) 
a value of Umax = 3 was chosen, the spectra can be eas¬ 
ily compared by normalizing the four data sets. For the 
quantitative comparison of spectra obtained for different 
driving amplitudes and, hence, different effective tunnel¬ 
ing amplitudes as in Fig. 3, however, normalization of 
each single spectrum is not a valid approach. For this 
reason, a large number of n^ax = 10^ pixels has been in¬ 
cluded in the calculation of the maximum optical density 
in Fig. 3(a). Here, the obtained data yields comparable 
results even for the zero-crossing of the effective tunnel¬ 
ing at A « 2.4. 

Despite the stark difference between the number of in¬ 
cluded pixels in the spectra of Figs. 2 and 3, their qual¬ 
itative behavior remains unchanged as shown in Fig. S3. 
Here, maximum optical density spectra obtained at a 
driving amplitude of K = 2.4 and lattice depths of (a) 
Vb = ll^-rec and (b) Vq = 13F^rec [compare Fig. 3(b) in 
the main text] are plotted over a wide range of values for 
^max- While, according to Eq. (SI), the absolute values 
of the maximum optical density necessarily decrease with 
increasing Umax it is evident that all important features in 
the spectra remain present without further data process¬ 
ing. This is the case even though for the total number 
of Up ~ 8.5x10^ analyzed pixels in each TOF image a 
value of nmax = 10^ pixels corresponds to the inclusion 
almost 12% of the absorption image in the calculation of 
the maximum optical density. 












